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I. INTRODUCTION 



Hadron scattering is presently one of the most intriguing process in high energy particle physics. Unlike the 
unification scheme which characterizes the electroweak sector of the standard model, some topical aspects of quantum 
chromodynamics (QCD) remain yet unknown and this has been a great challenge for decades. One point concerns 
some subtleties emerging from its running coupling constant, which entails that high energy hadronic phenomena 
have been classified into two wide and nearly incongrous areas, namely, large pr or hard processes and low px or soft 
hadronic physics. From a purely theoretical point of view (QCD), these phenomena are treated through perturbative 
and non-perturbative approaches respectively, and this renders difficult an unified formalism able to describe the 
totality of experimental data available on high energy hadronic interactions. The reason is that, despite the successes 
of perturbative QCD in the description of the hard (inelastic) hadronic scattering and also the successes of non- 
perturbative QCD in treating static properties of the hadronic systems, the scattering states in the soft (long range) 
region yet remains without a pure QCD explanation: Perturbative approaches do not apply and pure non-perturbative 
formalisms arc not yet able to predict the bulk of the scattering states. 

This soft hadronic physics is associated with the elastic and diffractive scattering, characterized, experimentally, 
by several physical quantities in both the elastic and inelastic channels, such as elastic differential cross section, total 
cross section, charged multiplicity distributions, average multiplicities and others JlJ. In spite of the long efforts to 
describe these data through a pure microscopic theory (QCD), our knowledge is still largely phenomcnological and 
also based on a wide class of models and some distinct theoretical concepts, such as Pomeron, Odderon, impact 
parameter picture, parton and dual models, Monte Carlo approaches and so on. 

However, presently, this phenomenological treatment of the soft hadronic physics plays an important role as a source 
of new theoretical insights and as a strategy in the search for adequate calculational schemes in QCD. The multiple 
facets associated with this phenomenological scenario have been extensively discussed in the literature and Ref. (|] 
presents a detailed outlook of the progresses and present status of the area. 

In addition to the intrinsic importance of high energy diffractive physics associated with our limited theoretical 
understanding, a renewed interest in the subject may be seen in the last years. This, in part, is due to the Hera 
and Tevatron programs, but also to the advent of the next accelerator generation, the Relativistic Heavy Ion Collider 
(RHIC) and the CERN Large Hadron Collider (LHC). In fact, with these new machines it shall be possible to 
investigate pp collisions at center-of-mass energies never reached before in accelerators, allowing comparative studies 
between pp and pp scattering at the highest energies, including both hard and soft processes. 

Presently, at this "pre-new-era" stage and due to the lack of an widely accepted unified theoretical treatment of both 
elastic and inelastic channels, it may be important to re-investigate ways of connecting these channels, looking for 
new information. Even if the treatment is essentially phenomenological, as explained before, the predictions shall be 
checked and may contribute with future theoretical (QCD) developments. 

To this end, in this work we shall investigate some aspects of both elastic and inelastic pp and pp scattering in 
the context of a particular phenomcnological approach. Our goal is to obtain simultaneous descriptions of some 
experimental data from both channels, that is, our primary interest concerns connections between elastic and inelastic 
hadron scattering. Accordingly we shall base our study on one of the most important principle of quantum field 
theory: Unitarity. 

For reasons that will be explained in detail in what follows, our framework is the impact parameter formalism 
(geometrical approach). At first, under geometrical considerations, we shall not refer to quarks/gluons or partons, 
but treat hadron-hadron interactions as collisions between composite objects made up by elementary parts, which we 
shall generically refer as "constituents" . At the end we discuss some possible connections between our results and 
the framework of QCD. Also, as shall be explained, our starting point is the description of physical quantities that 
characterize the elastic channel in pp and pp scattering. We then proceed to consider the inelastic channel through 
unitarity arguments and in a geometrical approach. 

Following other authors dfl 7 we shall express the "complex" (overall) hadron-p multiplicity distributions (inelastic 
channel) in terms of an "elementary" distribution (associated with an elementary process taking place at given impact 
parameter) and the inelastic overlap function, which is constructed from descriptions of the elastic channel data. The 
novel aspects concern: (a) quantitative correlation between the violations of the Koba-Nielsen-Olesen (KNO) scaling 
|^| (inelastic channel) and geometrical scaling Q (elastic channel) ; (b) introduction of novel parametrizations for the 
elementary quantities based on geometrical arguments and taking suitable account of the most recent data on contact 
interactions. With this general formalism, in addition to the description of the elastic data (even at large momentum 
transfers) , the hadronic multiplicity distributions may be evaluated and an excellent reproduction of the experimental 
data on pp and pp inelastic multiplicities is achieved |Q. We also present predictions at the LHC energies. 

The paper is organized as follows. In Sec. II we discuss the underlying phenomenological ideas, the data to be 
investigated and the strategy assumed. In Sec. Ill we present the theoretical framework connecting elastic and 
inelastic channels, the inputs from elastic scattering data, the novel parametrizations for the elementary processes, 
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the predictions for the hadronic multiplicities distributions and comparison with the experimental data. In Sec. IV 
we discuss in some detail all the results obtained and the physical and geometrical interpretations. The conclusions 
and some final remarks are the content of Sec. V. 

In what follows we shall represent the main physical quantities associated with hadron-hadron scattering (com- 
plex/overall system) by capital letters and those associated with constituent-constituent interactions (elementary 
process) by lower case. 

II. EXPERIMENTAL DATA AND PHENOMENOLOGICAL CONTEXT 

The broad classification in hard, soft (and also semi-hard) processes is based on the momentum transferred in 
the collision. On the other hand, depending on the physical process involved, high energy hadron scattering may 
also be classified into elastic and inelastic processes and the later, in diffractive (single and double dissociation) and 
non diffractive. Concerning both elastic and inelastic channels, one of the striking features that emerged from early 
experiments was the violation of the scaling laws, namely, the geometrical scaling in elastic scattering || and the 
KNO scaling in the inelastic events [^) . For this reason, our main interest in this work is to correlate quantitatively the 
above scaling violations and to discuss its phenomenological and dynamical aspects. To this end, before we present 
the underlying formalism and results, we discuss in this section the physical observables to be investigated and the 
reasons for our choices concerning phenomenology and strategies. 

A. Elastic channel 

The differential cross section is the most important physical observable in the elastic channel, since from it other 
quantities may be obtained, in particular, the integrated clastic cross section, o~ e i, and the total cross section, a tot 
(optical theorem) . The violation of the geometrical scaling may be characterized by the increase of the ratio o~ e i /o~tot 
with the energy at the CERN Intersecting Storage Ring (ISR) and at the CERN Super Proton Synchrotron (SppS) 
regions. 

The differential cross section yields the elastic hadronic amplitude, F (q, s), by 

do 



^=^(ff,-)l a (1) 
and this amplitude may be expressed in terms of the elastic profile function, T(b, s), by 

F(q,s)=ifbdbJ (qb)T(b,s), (2) 

where b is the impact parameterand ^Js the center-of-mass energy. 

As commented before, despite the bulk of models able to reproduce the differential cross section data at the 
ISR, SppS and Tevatron energies Q, an approach based exclusively in QCD is still missing. Obviously due to its 
soft character, a QCD treatment of the elastic scattering should be non-perturbative. Along this line, despite the 
difficulties, important results have recently been reached through the works by Landshoff, Nachtmann, Simonov, 
Dosch, Ferreira and Kramer (|||. The approach, based on the functional integral representation (QCD) and eikonal 
approximation, allows to extract a quark-quark profile function 7(6) (impact parameter space) from the gluon gauge- 
invariant two-point correlation function, determined, for example, from lattice QCD flio| |. Through the Fourier 
transform (analogous to Eq. (2) at the elementary level), the quark-quark scattering amplitude, f(q,s), may be 
obtained: 

f(q,s)=i [ bdbJ (qbh(b,s). (3) 



One possible connection with the hadronic scattering amplitude, Eq. (2), is by means of the Stocastic Vacuum 
Model (SVM) |J and some important results have recently been obtained jtl]]. However, presently, this theoretical 
framework still depends on some phenomenological inputs. Also, it is able to reproduce only the experimental data in 
the forward region and/or very small values of the momentum transfer and does not distinguish pp and pp scattering 
(dip region) , even at ISR energies jll) . 

Another way to obtain the hadronic amplitude from the elementary one is through the Glauber's multiple diffraction 
theory (MDT) []l2| and this plays a central role in our choices concerning phenomenology and strategies as discussed 
in what follows. 

Originally the MDT was applied to hadron-nucleus and nucleus-nucleus collisions |l^] and after to hadron-hadron 
scattering ]13[ |. The topical point which interest us is the allowed general connection between complex quantity 
(composite object) and elementary quantity (constituents). In the case of hadron-hadron collisions, the connection 
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between the hadronic amplitude (composite object) and the elementary amplitude (constituents) is also established 
through the eikonal approximation. In this approach the hadronic profile function, Eq. (2), is expressed by 



T(b,s) = l-e^ 6 ' 8 ) 



(4) 



where 



X (b,s)=C / qdqJo(qb)G A G B f 



(5) 



is the eikonal function, G a,b the hadronic form factors, / the elementaryf constituent - constituent) amplitude and 
C does not depend on the transferred momentum. The above notation ||l4f shall be useful for the discussion we are 
interested in. 

In spite of their simplicity, Eqs. (2), (4) and (5) are extremely useful. Recently, with suitable parametrizations for 
the form factors and with the elementary amplitude (quark-quark) extracted from a parametrization for the gluonic 
correlator through the functional approach (non-perturbative QCD), Grandel and Weise obtained good descriptions of 
the differential cross secion data for pp and pp elastic scattering at the IS R and SppS energies, but only in the region 
of small momentum transfer Jlq ). On the other hand, excellent descriptions of experimental data, including also large 
momentum transfers, have been obtained in a rather phenomenological context, through suitable parametrizations 
for G a-, Gb and / |16} |. Moreover, elementary amplitudes obtained through the SVM and the gluonic correlator from 
lattice QCD have been investigated and also comparisons with empirical analysis and model predictions have been 
discussed 

We understand that all these facts indicate that the impact parameter formalism (and the eikonal approximation), 
connecting the complex (overall) amplitude with the elementary amplitude (constituent-constituent), Eqs. (2), (4) 
and (5), is a very fruitful and simple approach in the investigation of the elastic hadron scattering. As shown, it seems 
also to be an adequate bridge between phenomenology and non-perturbative QCD. These conclusions constitute one 
of the foundations of our approach and, as discussed in what follows, the extensions to the inelastic channel shall be 
based on the general idea of connections between overall and elementary quantities in an impact parameter picture. 



Concerning scaling in the inelastic channel, the quantity of interest is the hadronic charged particle multiplicity 
distribution Pjv, normalized in terms of the KNO variable, Z — N(s)/ < N > (s), as 



The broader distribution observed at the SppS characterizes the violation of the KNO scaling, namely, $ = &(Z, s). 

As in the case of clastic differential cross section data, a wide class of models describes this behavior, as for example, 
dual parton pq |, fireball | p9| , two-component models p0| , and others. Also, hadronic processes have been extensively 
treated through Monte Carlo event generators and the Lund parton approach pl[-p3||. However, we observe that, 
despite some QCD inspired approaches and good descriptions of some soft processes, all these formalims and models 
are concerned exclusively with the inelastic channel and this is the topical point that distinguishes our strategy, as 
discussed in what follows. We shall also return to this subject in Sees. IV and V. 



Connections between Geometrical and KNO scalings were established a long time ago, by Dias de Deus ||, Lam and 
Yeung [Q. However, we are interested here in their violations and the central point is: Does one need a new connection 
when the two phenomena are violated at the SppS or can the two effects be correlated both phenomenologically and 
dynamically? We will argue that the latter alternative seems to be prevail. Specifically, our goal is to correlate 
quantitatively both violations in an analytical way and we shall show that, beginning with a formalism that describes 
quite well the violation of the Geometrical scaling (elastic channel input), it is possible to extend it and to describe, 
quantitatively, the violation of the KNO scaling in an analytical way. We stress that this strategy distinguishes 
our approach from all the other model and theoretical descriptions of elastic or inelastic scattering that treat these 
interactions separately, in an independent way or in distinct contexts. 

Since the connection between elastic channel — ► inelastic channel is our primary interest, the approach shall be 
based in direct analogy with the ideas discussed in Sec. II. A, that is, we consider hadron-hadron collisions as collisions 
between complex objects, each one composed by a number of more elementary ones. As an extension of the Glauber 
multiple diffraction theory, which connects hadronic and elementary elastic amplitudes, we shall consider the impact 
parameter formalism and also express the hadronic multiplicity distribution (complex system) in terms of elementary 
multiplicity distributions (constituents). The point is to describe a "wide" distribution (hadronic) by superimposing a 



B. Inelastic channel 



< N > (s)P N (Z) = 



(0) 



C. Strategies 
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number of narrower ones (elementary) p| . What we shall do here is to infer what these elementary distributions should 
be, in order to reproduce the experimental data on hadronic distributions and in the context of the impact parameter 
picture. He hope that, as in the elastic case, these information can contribute to further theoretical developments. 

III. UNITARITY AND IMPACT PARAMETER PICTURE 

Unitarity is one of the most important principles in quantum field theory. In the geometrical picture, unitarity 
correlates the elastic scattering amplitude in the impact parameter b space, T(b, s), Eq. (2), with the inelastic overlap 
function, Gi n (b 7 s), by 

2ReT(b,s) = |r(6,s)| 2 +G m (6,s) (7) 

which is term-by-term equivalent to j24j 

G tot (b,s) = G el (b,s) + G in (b,s). (8) 

For a purely imaginary elastic amplitude in momentum transfer space the profile function T(b, s) is real and in the 
eikonal approximation is expressed by 

T(b,s) = l-exp[-n(b,s)], (9) 
where f2(&, s) = Imx(b,s) in Eq. (4). With this, 

G m (b, s) = l- exp[-2 Q(b, a)] = a in (b, s) (10) 
is the probability for an inelastic event to take place at b and s and 

a in (s) = J d 2 b G m (b,s). (11) 

In this picture the topological cross section for producting an even number N of charged particles at CM energy 
yfs is given by 



0n(s) — J d 2 b <r N (b,s) = J d 2 b <j in (b,s) 



a N (b, s) 



(12) 



_C7j n (6, s 

where the quantity in brackets can be interpreted as the probability of producing N particles at impact parameter b. 

A. Hadronic and elementary multiplicity distributions 

We now introduce the multiplicity distributions for both an overall and an elementary processes in terms of corre- 
sponding KNO variables || and also the formal connection between these distributions. 

Representing the hadronic (overall) multiplicity distribution by $ and the corresponding KNO variable by 

z m as) 

Z -<N>(s) (13) 
where < N > (s) is the average hadronic multiplicity at we have in general 

$=<7V> (s)^44 =*{Z,a). (14) 

Now, let < n > (b, s) be the average number of particles produced at b and s, ip the elementary multiplicity 
distribution and 

- N(t) <«> 



<n> (b, s) 

a KNO variable associated with the elementary process taking place at b (and s). Then, in general, 

. . ^nQ>, s) 

tp =<n> (b,s) T7— r = (p{z,s). (16) 

Vin(0, S) 
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Both distributions are normalized by the usual conditions ||,|| 

p OO pOO 

I ${Z)dZ = 2 = / $(Z)ZdZ. (17) 
Jo Jo 

poo poo 

/ ip(z)dz = 2=1 (p(z)zdz, (18) 
Jo Jo 

The relationship between $ and ip then follows from Eqs. (10-12), (14) and (16): 

<N> (s) (tfh ^A y 
$ = ill <">(M)^ . . 

Jd?bG in (b,s) ■ [ J) 

Now, let us define a multiplicity function m(b, s) by the ratio 



„ < n > (b, s) , . 

-(M) = ^^, (20) 



so that Eq. (19) becomes 



«= J^bl(M (21) 

It is well known that connections between KNO and Geometrical scaling may be established if m(b, s) = m(b/R(s)) 
and also Gi n (b,s) = Gi n (b/ R(s)) ||, where R(s) is the "geometrical radius". In this case $(Z, s) is only a function 
ofZ. 

The general result (21) means that, once one has parametrizations for Gi n (b,s) and the elementary quantities ip 
(multiplicities distribution) and m(b, s) (multiplicity function) the overall hadronic multiplicity distribution may be 
evaluated. In this work we consider Gi n (b,s) from analyses of elastic pp and pp scattering data (taking account of 
geometrical scaling violation) and infer the elementary quantities based on geometric arguments, as explained in what 
follows. In so doing, we shall correlate quantitatively the violations of both KNO and Geometrical scaling in an 
analytical way. 

B. Elastic channel input: the BEL Gi n 

In the elastic channel, the breaking of Geometrical scaling is quite well described by the BEL behaviour, analytically 
expressed by the Short Range Expansion of the inelastic overlap function p4] , p5| 

G in (b, s) = P(s)exp{~b 2 /4B(s)}k(x, s), (22) 

with k being expanded in terms of a short-range variable x = b exp{—(eb) 2 /4i?(s)}, i.e. 

N 



k(x, s) = ^ 



n=0 



e exp{l/2} 



(23) 



The quantity in the bracket of (23) by itself exhibits GS for constant values of e 2 « 0.78, but k(x, s) doesn't because 
of the s-dependence of (>2n(s) and therefore Gi n (b, s) doesn't either (also because of P(s))- Each term in the bracket 
of (23) has a maximum value of 1 and the rapid convergence of the series reproduces data for all values of —t E (0, 14) 
GeV 2 with N=3. For k=l=P, we recover the Van Hove limit for cr e ; / cr tot =\- 4 ( 1 ^ ln2 ) ~ 0.1853 which is nearly attained 
at the ISR. The deviation of k from the constant value of 1, in particular the increase of with increasing s is 
responsible for the Edgier behaviour of Gi n (b, s), while increasing values of P(s) and B(s) make the proton Blacker 
and Larger respectively (BEL behavior of the inelastic overlap function) . 

Excellent agreement with experimental data on pp and pp elastic scattering is achieved ]2j| for the following 
parametrizations in terms of the Froissart-like variable y = ln 2 (s/so) with Sq = 100 GeV 2 

T) = 7 + ;y = 0.115 + 0.0009% ,24) 



G 



B(s) = 6.64 + O.OUy 
and 64 determined from 62 in some models |26| by 64 = <5|/4. 

C. Elementary hadronic process in a geometrical picture 

We now turn to the discussion of the elementary hadronic process, characterized by tp and m in Eq. (21). By 
construction, these quantities are associated to collisions of strongly interacting hadronic constituents. As commented 
before, due to the success of the geometrical models in the investigation of elastic hadron scattering (for example, 
the above BEL approach) and the presently lack of a pure QCD approach and/or Monte Carlo models to the subject 
(elastic/soft scattering states), we shall discuss what an elementary hadronic process could be in the geometrical 
framework and in an analytical way. Our arguments are as follows. 

In the geometrical approach an elementary process is a process occuring at a given impact parameter. Concerning 
contact interactions, experimental information is only available from lepton-lepton collisions, which is a process occur- 
ing in a unique angular momentum state and therefore also at a given impact parameter (zero in this case) . Although 
these processes can not be the same as collisions between hadrons constituents, it is reasonable, from the geometrical 
point of view, to think that some characteristics of both processes could be similar. The point is to find out or infer 
what they could be. 

For these reasons we shall consider the experimental data available on e + e~ collisions as a possible source of 
(limited) geometrical information concerning elementary hadronic interactions (at given impact parameter). We do 
not pretend to look for connections between e + e~ annihilations and pp and pp collisions but to extract from the 
former processes suitable information that allows the construction of the hadronic multiplicities (and the connections 
with the corresponding elastic amplitude) in an analytical way and in the geometrical context. This may be achieved 
for both tp and m in Eq. (21) through the following procedure. 

• Analytical relation between multiplicity function and eikonal 

First, in order to connect the multiplicity function m(b,s) and the eikonal Q(b, s) (and so Gi n (b,s) by Eq. 
(10)) in an analytical way, let us consider the very simple assumption that the average multiplicity at given impact 
parameter depends on the center-of-mass energy in the form of a general power law 

< n > (b fixed, s) oc E@ M . (25) 

We shall discuss this assumption in detail in Sec. IV. B.. 

Now, from Eqs. (10) and (11), exp{— 2 fi(6, s)} is the transmission coefficient, i.e. the probability of having no 
interaction at a given impact parameter, and therefore Q should be proportional to the thickness of the target, or as 
first approximation, to the energy Eqm that can be deposited at b for particle production at a given s. By Eq. (25) 
this implies 

< n > (b,s) oc ff(b,s). (26) 

Comparison of Eqs. (20) and (26) allows us to correlate the multiplicity function m(b, s) with the eikonal through a 
non- factorizing relation (in b and s): 

m(b,s)=Z(s)W(b,s), (27) 

with £(s) being determined by the normalization condition of the overall multiplicity distribution, Eq. (18). With 
this, Eq. (21) becomes 

r d 2 h G in (b,s) ( z \ 

Jd^bG m (b,s) [ °> 

where 

Jdb 2 G m (b,s) 

m ~ j dp G m (b,s)^(b, s y (29) 

Once the above analytical connection is assumed, the elementary hadronic process is now characterized by only two 
quantities, namely, the elementary distribution tp and the power coefficient 7. We proceed with the determination of 
these quantities through quantitative analyses of e + e~ data and under the following arguments. 
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• Elementary multiplicity distribution 

Because the elementary process occurs at a given impact parameter, its elementary structure suggests that it should 
scale in the KNO sense. Now, since experimental information on e + e~ multiplicity distributions shows agreement 
with this scaling J2?j, we shall base our parametrization for tp just on these data. In particular, it is sufficient to 
assume a gamma distribution (one free parameter), normalized according to Eq. (17), 

= 2^z K - 1 exp{-Kz}. (30) 

Fit to the most recent data, covering the interval 22.0 GeV < y/s < 161 GcV furnished K = 10.775 ± 0.064 

with x 2 /Ndf = 508/195 = 2.61 and the result is shown in Fig. 1. 

Concerning this fit, we verified that data at 29 and 56 GeV make the highest contributions in terms of \ 2 values. 
For example, if the former data are excluded we obtain K = 10.62 and x 2 I^df = 414/181 = 2.29 and if both sets 
are excluded then K = 10.88 and x 2 /N D f = 286/162 = 1.77. For comparison we recall that the DELPHI fit through 
a negative binomial distribution to data at only 91 GeV gives x 2 ' /Ndf = 80/34 = 2.35 (and x 2 /Ndf = 43/33 = 1.30 
through a modified negative binomial distribution) |27f |. However, we also verified that the above two values for K 
are not so sensitive in the final result concerning the hadronic multiplicity distribution, which is our goal (we shall 
return to this point in Sec. IV. A). For this reason and since we are only looking for experimental information that 
could represent contact interactions (geometrical point of view) we consider our first result shown in Fig. 1 as the 
representative one. 

• Power coefficient 

Finally, following Eq. (25), we consider fits to the e + e~ average multiplicity through the general power law 

<n> e+e -= A[^\\ (31) 

We collected experimental data at center-of-mass energies above resonances and thresholds and also the most recent 
data at the highest energies, covering the interval 5.1 GeV < yfs < 183 GeV |2^j2^]. Fitting to Eq. (31) yields A=2.09 
± 0.02, 

7 = 0.516 ±0.002 (32) 

with x 2 /Ndf = 409/46 = 8.89 and the result is shown in Fig. 2. We observe that this parametrization deviates 
from the data above y/s ~ 100 GeV and this contributes to the high x 2 value. However, as commented before, we do 
not expect that e + e~ annihilation exactly represent the collisions between hadrons constituents. The power law is a 
form which allows an analytical and simple connection between the multiplicity function and the eikonal as expressed 
by Eq. (27). In Sec. IV. B we present a detailed discussion concerning this power assumption and in Sec. IV. A and 
IV. C, we discuss the physical meaning of the differences between our parametrization and the experimental data on 
< n > e + e -. 

D. Results for the hadronic multiplicities distributions 

With the above results we are now able to predict the hadronic inelastic multiplicity distribution &(Z, s), Eqs. (28, 
29), without free parameters: Gi n (b,s) (and Q(b, s)) comes from analysis of the elastic scattering data (Eq. (10) and 
(22-24)); (p(z) and 7 from fits to e + e _ data, Eq. (30) and (31-32) respectively. 

We express $ in terms of the scaling variable Z = N / < N > where < N >= N(s) — Nq with Vo=0.9 leading 
charges removed. It is well known |3(| that such a subtraction improve the KNO curves for all measured data below 
the SppS Collider with the above value of Nq. This is completely equivalent to the Wroblewski |nj relation for the 
dispersion D = y/N 2 — < N > 2 =0.594 [< N > —Nq] with the same value of Nq. Values of iVo around 1 and the 
numerical value of the D vs N can be found by parton model arguments ]32] | for valence quark distributions. 

The predictions for pp scattering at ISR energies and pp at 546 GeV are shown in Figs. (3) and (4), respectively, 
together with the experimental data [ (33p4| ] . The theoretical curves present excellent agreement with all the data, 
showing a slow evolution with the energy at ISR for large Z and reproducing the KNO violations for large Z values 
at 546 GeV. In Fig. (4) is also shown the predictions at 14 TeV (LHC). 

IV DISCUSSION 

In the last section we obtained a quantitative correlation between the violations of both KNO and geometrical 
scaling. In the framework of the impact parameter picture, Sec. III. A, we only used four inputs, three parametrizations 
from fits to experimental data and one geometrical assumption, namely, 
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1. Gi n (b, s) from fits to elastic scattering data (BEL behavior); 

2. tp(z) from fit to e + e~ multiplicity distribution data; 

3. The geometrical assumption (27) concerning the multiplicity function m(b,s); 

4. The power coeficient 7, from fit to e + e~ average multiplicities data. 

In this section, we first investigate the sensitivity of each parametrization from fits to the experimental data in the 
output of interest, namely, the hadronic multiplicity distribution. Then, we discuss in some detail the geometrical 
assumption concerning the multiplicity function and the power law. Finally, based on this study, we outline the 
physical picture associated with all the results from both geometrical/phenomenological and QCD point of views. 

A. Sensitivity of the parametrizations 

First we observe that the power coefficient 7 in Eqs. (28-29) could, formally, be considered as a free parameter in 
a direct fit to the data on the hadronic multiplicity distributions and, in this case, it would not be necessary to take 
account of the e + e _ average multiplicity data. This, however, leds to a strong correlation between 7 and the other 
two inputs, Gi n (b, s) and ip(z). On the other hand, with our procedure, the values and behaviours of the three inputs, 
Gi„, f and 7, are rougly uncorrelated and this allows tests of the inputs by fixing two of them and changing the third. 
In what follows we perform this kind of analysis, beginning allways with the results obtained in the last section and 
considering, separately, a change in each one of the inputs. Since we are interested in the scaling violation we shall 
base this study in the results for the hadronic multiplicity distribution only at the collider energy. 



• Changing G in 

Among the wide class of models for Gi n [||, we shall consider a multiple diffraction model (MDM) and also the 
traditional approach by Chou and Yang, as a class of geometrical model (GM). The reason is based on the discussion 
in Sec. II. A. Also, as we shall show in Sec. IV. C, these models allows a suitable connection with the interpretations 
that can be inferred from our general approach. 

- Multiple diffraction model (MDM) 

This class of models is characterized by each particular choice of parametrizations for the physical quantities in 
Eq. (5), namely, form factors, Ga,b and elementary scattering amplitude / p4[ . In particular, Menon and Pimentel 
obtained a good description of the experimental data o pp and pp elastic scattering, above ^fs — 10 GeV through the 
following choices citemenonpimentel: 

G A (q) =G B (q) = J 1 (33) 

(! + &)(! + £) 

The parameters a 2 and (i 2 are fixed and the dependence on the energy is contained in the other two parameters: 

C(s) =&exp{UlHs)?}, (35) 



af s) =^Ms)} £2 , (36) 

where £j, i = 1,2,3,4 are real constants. The reason for these choices and physical interpretaions are extensively 
discussed in [ |36[ . With these perametrizations the opacity function, 

0(M) = Im X (M), (37) 

is analytically determined and then the inelastic overlap function through Eq. (10). 
- Geometrical model (GM) 

In the geometrical approach by Chou and Yang, the essential ingredient is the convolution of form factors in the 
impact parameter space |p7[ . However, in the context of the multiple diffraction theory, it can also be specified by 
the following choices [ |l4"| ]: 

G A (q) = G B (q) = l — , (38) 

2^[I + (^)] 2 
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/ = !• 



(39) 



In Ref. |38| the parameters /i and C were determined through fits to elastic pp data at y/s = 23.5 GeV and pp at 
546 GeV. Following the authors, we consider the parametrizations 

C(s) = a\ + a 2 lns, = b\ + 6 2 lns. (40) 

A* ( S J 

With the above double pole parametrization for the form factors, the opacity, Eq. (5), is analytically determined 
and so the inelastic overlap function, Eq. (10). 
- Results 

The results for the inelastic overlap function at 546 GeV are shown in Fig. 5, from both the MDM and the 
GM, together with the BEL Gi n for comparison. In what follows, we shall use the following notation for these 
parametrizations G^ DM , G^ M , and G? n EL , respectively. 

We then calculate the hadronic multiplicity distribution, Eqs. (28-29), at this energy, by fixing both 7 = 0.516, Eq. 
(32), and the gamma parametrization for the elementary multiplicity <p(z), Eq. (30), and using Gf^ DM and Gf^ 1 ■ 
The results are displayied in Fig. 6 togheter with that obtained with Gf EL (Fig. 4) for comparison. 

We observe that, for central collisions (small b), Gf EL and Gf^ 1 are very similar, but Gf^ DM has higher values 
(Fig. 5a). This leads to the differences in cf>(Z') at high multiplicities, as can be seen in Fig. 6. In the same 
way, the differences between Gf„ DM , G^ M , and Gf EL at large b (Fig. 5b) originate the differences in cf>(Z') at 
small multiplicities. In all the cases the physical picture is that large multiplicities (large Z ) occur for small impact 
parameters while grazing collisions (large b) lead to small multiplicities, as one would have naively expected. 

An important conclusion is that, with 7 and ip(z) fixed, the hadronic multiplicity distributions obtained with 
qMDM ^ qGM ^ anc } QBEL re p roc i uce the experimental data quite well. We shall return to this point in Sec. IV. C. 



• Changing the elementary distribution ip 

As a pedagogical exercise, we shall consider only an early parametrization introduced by Barshay and Yamaguchi 

. , 817T 2 3 . 97T 

Vby{z) = —z 3 exp{-—z 2 }. (41) 

This function was used in the analysis of e + e~ multiplicity distributions at lower energies and, as can be seen in Fig 
7, does not reproduce the data at higher energies as well as the gamma parametrization. 

As before, we now proceed by fixing both 7 = 0.516 and Gf n EL and using the above parametrization for the 
elementary distribution. The result for the hadronic multiplicity distribution at 546 GeV is shown in Fig. 8, together 
with the result obtained with the gamma parametrization for the elementary process (Fig. 4). The broader width of 
iPby(z) as compared with that of the gamma distribution, is directly reflected in the hadronic multiplicity. Despite the 
differences between the two parametrizations for the elementary process the final result for the hadronic distribution 
with if by can yet be considered as a resonable reproduction of the experimental data. 



• Changing the power coefficient 7 

Finally, we consider different parametrizations for the e + e~ average multiplicity data in the interval 5.1 < y/s < 183 
GeV ,but under the assumption of the power dependence. We shall discuss this assumption in the next section. 
First we consider the naive parametrization based on the thermodynamic model (see next section) 

<n> e+e -=2.20[V^]°- 500 . (42) 

For the above ensemble of data one obtains \ 2 /DOF = 209/48 = 4.35. 

Second, and more importantly, we shall investigate the effect of the data at the highest energies, which are not 
reproduced by our original parametrization, as can be seen in Fig. 2. To this end, we consider only the data above 
10 GeV (25 data points) and the general power law parametrization. With this procedure we obtained 

<n> e+e -=3.46[^°- 396±0 - 008 . (43) 

with x 2 /DOF = 27/23 = 1.7. 

The result is displayied in Fig. 9 together with Eq. (42) and our original parametrization, Eqs. (31-32). We observe 
that, concerning e + e~ average multiplicity Eq. (43) brings information from data at high energies (roughly above 50 
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GeV), while the original parametrization, Eqs. (31-32) is in agreement with data at smaller energies (below ~ 100 
GeV) and the same is true for the parametrization with Eq. (42). 

As before, we now calculate the corresponding hadronic multiplicity distribution by fixing both the gamma 
parametrization for the elementary distribution, Eq. (30), and the G^ Bi , Eqs. (22-24), and considering the three 
parametrizations for the average multiplicity, Eqs. (31-32), (42) and (43). The results at 546 GeV are shown in Fig. 
10. 

We conclude that, in the context of our approach with the fixed inputs Gf n EL and gamma parametrization for 
f(z), the information from the e + e~ average multiplicites at high energies with the power-law does not reproduce the 
hadronic multiplicity distribution. That is, the elementary average multiplicity distributions in hadronic interactions 
must deviates from the < n > e + e - (s) as the energy increases, roughly above ~ 50 — 100 GeV. We shall discuss the 
physical interpretations of this result in Sec. IV. C. 

B. The multiplicity function and the power assumption 

We now turn to the discussion of a crucial assumption in our approach, namely, that the elementary average 
multiplicity at fixed impact parameter collisions grows as a power of the center-of-mass energy. To this end we shall 
first briefly recall some aspects of the power-law in hadron-hadron and e + e~ collisions, both in experiment and theory, 
and after, based on these ideas, we shall present a discussion concerning the use of this assumption in our approach 
and also the meaning of the multiplicity function. 

From the early sixties cosmic ray results on extensive air showers, at energies Eiab < 10 6 - 10 7 GeV, led to empirical 
fits of the type < TV >cx cx [-v/s] 1 ^ 2 (see for a review). A general power law with the exponent as a free 

parameter was used a long time ago, in order to allows analytical connections in analysis of cosmic ray data fl4l| . 
Also, in the beginnig of sixties, these investigations introduced the concept of inelasticity E9). This comes from the 
observation that the energy effectively available for particle production could not be identified with the cm. energy, 
as believed before 1953 (Wataghin, Fermi, Landau), but only with a fraction of it: 

W = ky/s. (44) 

The remained (1 — k)y/s was associated with the early named "isobar" system, presently known as leading particle. 

From the theoretical side, the power dependence emerged in the context of statistical models (Fermi, Pomeranchuck) 
and hydrodynamical models (Heisenberg, Ladau) pOfl . For example, taking account of the inelasticity, in the Landau 
model, the fact that the averaged multiplicity is proportional to the total entropy leads to the result f42]j 

< N >oc fc 3 / 4 ^] 1 / 2 . (45) 

Dependences on s 1 / 2 is characteristic of the Heisenberg and Pomeranchuck models and even s 1 / 8 appears in the 
Landau model, when viscosity is taken into account pPf . In the context of termodynamic models, a universal formula 
was discovered for proton targets and for energies below ~ 50 GeV: Data including 7, w, N and p collisions with p 



were quite well reproduced by < N >= 1.75 s 1 / 4 |43|. Concerning e + e data on average multiplicity this model 
suggested < n >= 1.5 s 3 / 8 and pure fits to low energy data furnished < n >= (2.2 ± 0.1)s 25±O Q1 , and also 
< n >= (1.73±0.03)s a34±0 01 Moreover, the power-law, with the exponent 1/4, was successfully used in the 



context of the parton model, either connecting KNO and Bjorken scaling |46j] or treating directly the violation of the 
KNO scaling @. 

The power-law may also appear under more general arguments. For example, suppose that an intermediate state 
(fireball) of invariant mass M oc y/s decays into two systems each of invariant mass Mi = M/c, where c is a constant. 
Suppose also that similar processes continue through some steps (sucessive cluster production) until the masses reach 
a value Mo (some minimum ressonance mass). It is easy to show that the final multiplicity reads p2f 

n oc [ % /i] ln2 / lnc . (46) 

For example for c = 4 — > n cx s 1 / 4 . The exponent 7 (our notation) may be inferred from c = 2 1 / 7 , so that higher 7 
values imply in higher splited masses in each step (for 7 = 0.516 — > c ~ 3.8). 

Based on the above review, we see that the power-law is characteristic of several analysis of experimental data on 
hadron-hadron and e + e _ collisions and also several theoretical approaches and models. Now we shall discuss this law 
in the context of our approach. 

First let us stress that in our formalism this assumption concerns an elementary hadronic process taking place at 
fixed impact parameter b. Thus, it does not pretend to represent the average hadronic multiplicity < N(s) >. Also, 
we used e + e~ data only as a possible source of information on contact interactions (fixed b) and therefore the power 
assumption does not pretend to represent the average multiplicity in e + e~ collisions. This is a subtle point in our 
approach and we would like to discuss it in some detail. 
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The main reason for the power assumption was to obtain an analytical and simple connection between the multi- 
plicity function m(6, s) and the eikonal, Eq. (27), which allows the general analytical connection between the elastic 
and inelastic channels. Since it is typical of several kind of collisions, as reviewed above, it is not unreazonable that 
it could represent an elementary hadronic process taking place at fixed impact parameter. Just for its elementary 
character (at given 6), there seems also to be no reason to include any inelasticity effect (leading particle) in the basic 
assumption represented by Eq. (25). That is, it seems reasonable that < n > (6 fixed, s) may be just proportional 
only to ET CM . 

The multiplicity function m(6, s), as defined by Eq. (20), connects the hadronic and elementary (at given 6) average 
multiplicities. With the power assumption and the geometrical arguments of Sec. II. C, m(b, s) may be expressed in 
terms of the eikonal and the power coefficient 7. The subtle point in our approach is that, since by definition 771(6, s) 
is proportional also to the average elementary multiplicity at given 6, the coefficient 7 was determined by fit to data 
available on contact interactions. In this sense, the model "imposes" the power-law and the e + e~ data are supposed 
to provide the limited, but possible, information on contact interactions. 

These considerations may allow to infer a distinction between e + e~ average multiplicity and what this quantity 
could be in an elementary hadronic process. Specifically, we showed in Sec. IV. A that data on the average multiplicity 
in e + e~ collisions, presently available above 5 GeV , can not be reproduced by the power-law. For example, a second 
degree polynomial in Ins gives a quite good fit to all the data above 5 GeV: 

< n > e+e - (s) = 0.0434 + 0.7751ns + 0.1681n 2 s (47) 

with x 2 /DOF = 145/45 = 3.2. However, besides this parametrization does not allow the analytical connection 
with the eikonal, we showed that, with the power-law, the behavior of < n > e + e - (s) at energies above ~ 50 GeV 
does not lead to the description of the hadronic multiplicity distribution. In other words, in the context of our 
approach, the increase of the elementary average multiplicity with energy in hadronic collisions must be faster than 
that observed in e+e~ collisions. This is not the case at lower energies, since the power-law with 7 = 0.516 gives a 
satisfactory description of the e + e~ data. In the next section we discuss the physical interpretations associated with 
these observations. 

C. Physical picture 

Based on the results of Sees. II and III, we now discuss the physical picture associated with the scaling violations, 
specifically, with the evolution of the hadronic multiplicity distribution $(Z') from the ISR to the collider and LHC 
energies, Figs. 3 and 4. 

From Eq. (21) the hadronic multiplicity $ is constructed in terms of Gi n and the elementary quantities if and m. 
In our approach, if scales and so does not depend on the energy. The multiplicity function m(6, s) is connected with 
G m through Eqs. (10) and (27), 

m(6, s) = £(s){hi[l - G m (s, b)}}\ (48) 

where £ comes from the normalization condition (29). Both £(s) and m(b, s) depend on the power coefficient 7, which 
is a constant determined from the fit through Eq. (31). Therefore, the evolution of the hadronic distribution with 
energy comes directly from Gi n (b, s) and depends also on the value of the exponent 7. This exponent, in turns, comes 
from the elementary average multiplicity dependence with the energy, Eq. (31), and therefore is associated with the 
effective number of colliding constituents in the hadronic process. 

Based on the above observations, the physical picture that emerges is that the energy evolution of the hadronic 
multiplicity distribution is correctly reproduced by changing only the overlap function, without tampering with the 
underlying more elementary process (f) . The geometrical evolution of the constituents of the hadron is responsible 
for the energy dependence and not the dynamical interaction itself. This is what one would expect if the underlying 
interaction is unique (QCD) but the relative importance of the constituents involved in collisions changes with energy 
(indicated by the exponent 7). 

We showed that with the power assumption, the information from e + e~ data above, say, ~ 100 GeV leads to 
an understimation of the hadronic multiplicity ditribution (Fig. 10). This means that the average multiplicity in 
an elementary hadronic process must increase with energy faster than that associated with e + e~ collisions. This 
result seems quite reasonable since, in a QCD guided approach, we expect different contributions from gluons/quarks 
interactions than those associated with lepton-lepton collisions. As the average multiplicity increases, the relevance of 
the original parton decreases, so that at high energies e + e~ can serve as a good first guide to quark-quark, quark-gluon 
and gluon-gluon multiplicity distributions. In a parton model (following QCD), this effect above ~ 100 GeV may be 
interpreted as the unset of gluons interactions |32]]. The faster increase represented by our power-law with 7 = 0.516 
(Fig. 2) may be atributed to the fully development of the gluonic structure, rather than the quark (valence) structure. 
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These "microscopic" interpretations are also directly associated with the BEL behavior, since its origin may be 
traced either to gluon interactions in the eikonal formalism Jl7],[l8|], or to the increased size of spot scattering in the 
overlap function formalism [ fl9"|] . 

As commented before, a novel aspect of this work concerns the simultaneous treatment of both the elastic and 
inelastic channels. Specifically, we started from elastic channel descriptions (pp and pp differential cross sections) and 
extended the results to the inelastic channel (multiplicity distributions). In this sense, we expect that the physical 
picture from both channels should be the same. Besides the microscopic interpretation associated with the BEL 
Gin(b,s), even if we consider the naive models represented by the MDM and the GM, discussed in Sec. IV. A, the 
same scenario emerges. In fact, in both models the elementary interaction, represented by the elementary elastic 
amplitude f(q, s), does not depend on the energy, Eqs. (34) and (39). The energy dependence is associated with the 
form factor G(q,s) and the "absorption constant" C(s). The former, through the associated radius, 

i* 2 ( S ) = -6^| ?2=0) (49) 

describes the expansion effect (geometry). The latter is associated with absorption (blackning) in the context of the 
geometrical (Chou-Yang) model and to the relevant number on constituents in the context of the multiple diffraction 
theory [Q. herefore, we also conclude that the elementary interaction is unique (does not depend on the energy), but 
the geometrical evolution of the constituents and its relevent number in collisions changes as the energy increases. 

V. CONCLUSIONS AND FINAL REMARKS 

The underlying theory of the hadronic phenomena is QCD. As commented in Sees. I and II, depite all its successes, 
the theory has presently some limited efficiency in the treatment of soft hadronic processes, meanly related with 
unified descriptions of physical quantities from both elastic and inelastic channels. Moreover, some QCD approaches 
are based in extensive Monte Carlo calculations and concerning this point, we understand that, although these 
techniques represent a powerfull tool for experimentalists, it is questionable if they could really be the adequate and 
final scenario for a theoretical understand of the hadronic interactions, mainly if we think in connections with first 
principles of QCD. 

At this stage, it seems that phenomenology must play an important role to bridge the gap, or, at least, to indicate 
or suggest some suitable calculational schemes for further theoretical developments. On the other hand, all the 
phcnomcnological approaches presently available, have also very limited intervals of validity and efficiency in the 
treatment of hadronic processes at high energies. One of the serious limitations of the geometrical approach is the 
difficulty to directly connect its relative efficiency with the well established microscopic ideas (QCD). However, it has 
not been proved that this direct connection can not be obtained. 

In this work, making use of the unitarity principle and in the context of a geometrical picture, we obtained 
analytical connections between physical quantities from both clastic and inelastic channels. In particular we correlated 
quantitatively the violations of the geometrical and KNO scalings in an analytical way. The physical picture that 
emerges from both channels, for pp and pp collisions above ~ 10 GeV is the following. The dependence of the physical 
quantities with the energy (elastic differential cross section and inelastic multiplicity distributions) is associated with 
the geometrical evolution of the constituents and the relative importance of the constituents involved in the collisions. 
The underlying elementary process or interaction does not change with the energy. This is in agreement with what 
could be expected from QCD. 

With this kind of approach the correct information extracted from the elastic channel is fundamental. Our prediction 
at LHC energies was based on extrapolations from analysis at lower energies and so has a limited character. This 
observation, and obviously other considerations regarding different models, point out to the importance of complete 
measurements of physical quantities associated with the elastic channel at the LHC, that is not only total cross 
sections but also the p parameter and differential cross sections at large momentum transfer. 

Based on the limitations referred to in this section, we do not pretend that the forms we inferred for the hadronic 
constituent-constituent collisions, multiplicity distribution and average multiplicity, are a conclusive solution. How- 
ever, we hope that, at least, they can bring new information on what some aspects of an elementary hadronic process 
could be. 
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Figure Captions 



Fig. 1. The KNO charged multiplicity distribution for e + e annihilation data p7| , [28| and the fitted gamma 
distribution, Eq. (30) (dashed). 

Fig. 2. The average charged multiplicity for e + e~ annihilation data p8j , p9| and the fitted power law, Eq. (31). 

Fig. 3. Scaled multiplicity distribution for inelastic pp data |33| at ISR energies compared to theoretical expectations 
using Eqs. (28-29). 

Fig. 4. Scaled multiplicity distribution for inelastic pp data |}4j at 546 GeV compared to theoretical expectations 
using Eqs. (28-29) (solid) and predictions at 14 TeV (dashed). 

Fig. 5. Inelastic overlap functions for pp collisions at 546 GeV, predicted by the multiple diffraction model (MDM), 
geometrical model (GM) and the short-range-expansion, black-edge-large approach (BEL): (a) central region; (b) 
large distances. 

Fig. 6. Same as Fig. 4 with the three different inputs for the inelastic overlap function. Same legend of Fig. 5. 
Fig. 7. Same as Fig. 1 with the Barshay-Yamaguchi parametrization, Eq. (41). 

Fig. 8. Same as Fig. 4 with two different inputs for the elementary multiplicity distributions: gamma function, 
Eq. (30) (solid) and Barshay-Yamaguch parametrization, Eq. (41) (dot-dashed). 

Fig. 9. Same as Fig. 2 with three power-law parametrizations: 7 = 0.516 (solid), 7 = 0.500 (dotted) and 7 = 0.396 
(dot-dashed). In the last case only data above 10 GeV were fitted. 

Fig. 10. Same as Fig. 4 using the three different parametrizations for the elementary average multiplicity (Fig. 9 
and same legend). 
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